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INTRODUCTION 


A computational  propram  Inis  boon  prepared  by  tin*  Information 
Processing  Division  to  compute  a particular  Roradiat  ion  Funot ion  w(xf), 
Trans  for  Function  W(w,P),  Impulse  Response  w(t,pl,  and  Reradiatod 
Waveform  r.(t,P)  ns  ilo l'inod  below  In  terms  of  the  incident  plane  wave 
pulse.  This  IBM  04  Program  designated  USI,  Program  No.  0*'3U , is  in 
Fortran  tl  language  and  is  described  in  Appendixes  A and  R.  Similar 
computational  programs  are  described  in  USE  'IVehnieal  MemoraV.dum  No. 

2242-156-67  and  2242-157-67. 

t 

THEORY 

Reference  (a)  contains  a description  of  the  mathematical  model 
construeted  nnd  t.lif  theory  behind  cons  Idering  rc't’leetion  as  a re- 
rad inf  1 o.i  phenomenon . 

!' V ; program  computes  for  integer  values  of  \) 

Ai\  i '» 

,a)  nx  sa  w (x)  N ■„  r 

<b)  ''  (•■ »f'  ' ' P crTl 


^ i 
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(c)  kw  (t ,Y> 

(d)  V,  (t.P) 


who  c 


(O  w(X) 
(2)  w(w,f) 


W - 

S3 

(3)  w(t,f)  - 

Hen(x)  -= 


(-0° 

a,az 

n*  dx* 

Laguerre  polynomial 

I 


O'  4 vTtt 

di"  v / 


“ Hermite  polynomial 
■oG 


(M  e(t,P)  - | 

f(t) 


-«o 


( AU|  - o f |t|  r 7 


COMPUTER  PROGRAM  DESCRIPTION 

A nomenclature  listing  for  USL  Program  No.  0836  is  Appendix  A, 
the  - ow  chart  is  Appendix  D,  and  the  IBM  70-'*  Fortran  II  Program  i; 
Appendix  C. 


The  basic  input  data  dock  required  by  the  program  consists  f 
cards . 


feu  1 
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1-8 
y-16 
17-24 
25-32 
33-40 
4 1-48 


“1 

*2 

X1 

x2 

c 

V 

0 

ISKP  (set  equal  to  zero  to  compute 
Reradiation  Function) 

JSKP  (set  to  zero  to  compute  Transfer 
Function) 

KSKP  (set  equal  to  zero  to  compute 
Impulse  Response  & Reradiated  Waveform) 
For  long  jobs  requiring  the  use  of  a 
dump  tape  at  least  one  of  the  above 
option  variables  should  not  be  set 
equal  to  zero. 

Initial  value  of  u) 

Maximum  value  for 

Initial  value  of  t 

Initial  value  of  k 

KK  (if  set  equal  to  zero,  initial 

value  of  k will  be  computed) 

Maximum  value  of  k 


A t 


Increment  of  u) 


A k 


Components  of  X 
Components  of  N 
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Table  l (coat’d) 


Curd  Ko.  Cols. 


Contents 


4-11 

12-19 

20-27 

?6-35 

36-43 

44-51 

52-59 


Number  of  values  of  Reradiation 
Function  to  compute 

AT 

Maximum  value  for  t 

u>o 

A u) 

Maximum  value  of  X 
Initial  value  of  X 

6) 


FORMATS 

Card  No.  1 - Fo  mat  6F8.3,  413 

2 - Format  4pfi.3,2X,I2,2X,F6.3 

3 - Format  9F$ . 3 

4 - Format  13,  7F8 . 3 


Tape  Units  Required 


Tape  Unit  No. 


Tape  Identification 


6 

7 


Data  input 

Values  for  Re radiation  Function, 
Transfer  Function,  and  Impulse 
Response 

Calcomp  Plotter  containing 
values  for  Reradiation  Function 
Reradiated  waveform  array  (k,t) 
Transfer  Function  Array  (k 


US  I.  Tech.  Memo. 
2242-111-67 


Tape  Units  Required  (eont'd) 
Tape  Unit  !'o. 
fi 

( SS-5  must  be  down  to  dump 
No  other  sense  switches  are 


Tape  Identification 
Impulse  Response  Array  (k,t) 
Dump  Tape 

used 


Subroutine  Amp  required  - computes  values  of  A array  referred  to  unde 
equation  h. 


PROGRAM  OOTPUT 


Tape  #4  cental  ns  : 

(1)  The  values  for  the  i arruy  plus  the  corresponding  values 
for  the  Reradiation  Function  according  to  Format  (IX,  F10.5,  bX,  F10.S) 

(2)  The  values  for  the  produet  of  k and  w plus  the  corresponding 
values  of  the  Transfer  Function  according  to  Format  (IX,  F10.5,  5 X, 
F10.5) 

(3)  The  va  ues  for  t/k  plus  the  corresponding  values  for  the 
Impulse  Response  according  to  Format  ( IX,  FlO.5,  5X,  F10.5) 

Tape  tr*  contains: 

The  va  ues  for  the  Rcradiat.ion  F’unction  (Calcomp  Plotter  tape) 

Tape  /¥>  contains: 

The  Reradiatcd  Wavefonn  Array  (k,t)  according  to  Format  (FIO.1^) 
Tape  ,7  7 contains: 

The  Transfer  Funoti  n Array  (k,w)  with  Format  (FlO.b) 

Tape  ,'fS  contains: 

The  Impulse  Response  Array  according  to  Format  (FlO.b) 
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T&petfU  ; r-  a runp  tape . 

Notes:  This  projrHn  oonta  ns  options  to  compute  or  not  to  compute 
any  of  the  fu  >ct : one  mentioned  ab^ve.  Tapes  Unit  Nos.  6,  7,  and  8 
ran  be  used  jo.  .. nput  to  USL  Program  No.  O0OU , "Hepresontation  of 
Surfaces;  A C‘Tiput.er  Program  to  Plot  Contours  and  Draw  Perspective 
Views",  by  '•i.lwnri  hoards worth,  Jr. 


, SlfMMARY 

* An  IBM  704  Fortran  program,  USL  Program  No.  0838,  has  been 

written  to  eonputo  a particular  Reradiation  Function,  Transfer 
Function , Impulse  Response,  and  Rerad  ated  Waveform  in  terms  of  the 
incident  plane  wave  pulse. 


0.0. 

D.  A. 
Mathematician 


ri 

4 
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APPENDIX  A 

NOME  CLATURE  LISTING  FOR  USL  PROGRAM  NO . 0838 


S(I) 

k1  *-k> 

Z(I) 

\TsuT 

RERAD  (I) 

Element,  of  Rerad  ation  Function  Array 

TRFER  (I) 

Element  of  Transfer  Function  Array 

AXW  ( LM , I ) 

k.w 

RSSP  ( LI-1, 1) 

Element  of  Impulse  Response  Array 

RATIO  (LM,l) 

t'k 

GSUM  ( LM, J ) 

Element  of  Re radiated  Waveform  Array 

AT 

*1 

AO 

Ro 

c 

X' 

xi 

X? 

*2 

c 

c 

V 

V 

N 

w 

w 

WMAX 

Maximum  value  for  w 

T 

Initial  va  ue  for  t 

A K 

k 

AKMAX 

Maximum  value  for  k 
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r 


31 

B2 


B4 

3? 

56 

37 


B8 

310 

Bl? 

One^a 

DELTA 

TT 

TAV 

PHI 


III 

TMAX 


components  of 


components  of  N 


A x2 

A 

Increment  of  w 


A 't 

A k 

A). 


A w 

Maximum  value  for  'T 
Initial  va  ue  of  ^ 

4 


Number  of  values  of  Reradiation 
Function  to  compute 

Maximum  value  for  T 


9 


■ 

f* 
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APPENDIX  P 

FLOW  CHART  FOR  USL  PROGRAM  NO.  O838 
"RERAOIATIO::  FUNCTION,  TRANSFER  FACTION,  L’^IPULSE  RESPONSE  (CASE  2)" 


c 

c 


r 


MfWAriATICN’  FUNCT  ION. TRANSFER  H.NC  T JON  , \ mPUL  SE  R{ 5PCNSF  (CASF  2> 
r.  A.STPF^'SKY 

rivFNMCN  / ( 500  ) »RFP AC  ( > »TRF  ER< 50.40 ) «AHW( 50 ,50) .RATIO <40,50 ) #R 

1FSP(5p»50)  ,SU'M<5P)  ,PSUM(40)  .5(400)  ,GSUM(40,50)  ,R<1000)  ,I0UMP<1A> 

C I M£  N4 I 0^  HUF  F f R < 1 024 ) ,XAXlS(500)»YAXlS<400) 
stcjn  alf  *oh3h 

STECC  A|_f  • 

WRITF  OUTPUT  T APE  4,9501 

<>501  format  f i m i t 

REAP  INPUT  TAPF  3. 9502, Cl 

9502  FORMA T ( A5 ) 

IF  (TOjN-t'I  ) 940 1.9504,940* 

9503  PALSF  6 

9504  WRITE  OUTPUT  TAPF  4.9502.CI 

WRITE  OUTPlT  TAPF  4.9505 

9505  FORMAT  ( 10X32F'C.A.STRFN'S<Y«PCCM  3126,00^  2242) 

rfac  IFPUT  TAPF  3,100.A1,A2,X1.X?,C,V»N.ISKP.JSKP.KSKP 

100  FORMAT  (6FP.3.4I  3) 

READ  INPUT  TAPF  *,101 .W.WMAX.T.AK.KK.AXMAX 

101  FORMAT (4FH.3.2X,I?,2X.F8.3) 

READ  INPUT  TAPf  *.10‘'.P1»P2.P1,P4,P5,P6.P7»R«»R12 

103  format (9F8. 3) 

RFAC  INPUT  TAPE  3.10A.N1.P10.TMAx,AMFGA,DFLTA,TT.TAU»PM 

104  FORMAT ( I 3.TF». 3) 

w i »w 

A< 1 ■ AX 
T 1 -T 

tal i«tau 

PlF»3. 1415 
PCP=2.50O6 
CFG-1AP./PIE 
NP 1 *N«  1 
N P 2 « N ♦ 2 

Su* (1 ) *0. 

5lm (2) * 1 • 

IF (ISxP) 195.102.194 

102  J«1 

105  5(J)«<X1/A1)44?.(X2/A2)**2 
2X»S(J) 

Z fJ)«sORTF (ZX) 

IF (N) 146,156, 157 

156  S I G \ ■ 1 » 

GC  TO  158 

157  SIGN* <-l ) «*N 

1 58  RFRAP ( J) «SIGN»2.0*PIF«A1»A2 
IF  (Nl-J) 546,546, 1 75 
175  J*J*l 

X 1 «X 1 *F  5 
X2*X2*P6 
CC  TO  lO’' 

546  WRITE  OUTPUT  TAPE  4,551 

551  FORMAT ( 1 X36MZ  REmAHJATION  FUNCTION) 

WRlTt  OUTPUT  TAPI  4.552»(Z(I) «RERAO(H  ,I«1,N1) 

552  FORMAT  1 1X.F 10.n,5X,F1C,5) 

WRITE  OUTPUT  TAP)  4,553 

553  Format (///) 

194  IMsl 


I if 


IF  (jSxP)2B6*20b»2F6 
200  Jcl 

548  MJITt  OUTPUT  TAPE  4*554 

554  Fcbmt<IX34H'u  TRANSFER  FUNCTION) 

IF  IK.O  ?05«202*?05 
202  P 1 *P  1 -C*P3/V 
P2»B2-C*t'4/V 

CP*(A1*PI)*»2MA?«P2)**2 
AKsscDTF  < CH > /C 
205  lF(AK-AKMAX)207, 207*285 
207  IF (W-*MAX)210.210.250 
210  AkmLv,j)  *AK«H 
R?  = AK’* <LN» J) **2 
R3  = -R 2/2. 

CCEF  UFXPF  (R3) 

IF (N)? 30 *2 30 *215 
215  CC  225  I *3 »NP2 
AJV3»I-3 
I M 1 « I - 1 


lN2«I-2 

PAR«2.*AIN3*l.-k? 

PART»AIM3««2 

SU^ (I > «-PAR«SUV ( IM1)-PART*SUV (IN?) 

225  CONTINUE 
230  NPRCD-1 

IF(N)?P0*280*260 
260  CC  275  La  1 *N 
NE^L«L 

NPRCO»NPRCD«NEv*L 
275  CONTINUE 

280  PRCDoNPROD 
lF(N)281*281*282 

281  SIGNbI. 

GO  TC  283 

282  5IGN»(-1>*«N 

283  TRFER  <LM*J) cCCFF1*5UM(NP2)«SIGN/proo 
J»  J*  1 

V*sV**Pfl 
CC  TC  207 
250  N2«J-1 
N4  «L  V 

WPITt  OUTPUT  T APf,  4.555*  ( (AKV*(LN,I)  *tRFER<LN*I)  »I«1*N2)  *LM«N4.NA) 
555  FCRNATI1X.F10.5.5X.F10.5) 

533  V*R  I TE  OUTPUT  T APf  7 . 5 1 5 . ( ( TRFER  (LN  , I ) * I ■ 1 *N2 ) *1  M»N4  *N4 ) 

515  FCRNAT (F10.5) 

J«1 

W iW  1 
L v s L m ♦ i 

Ak=AK*P12 
GO  TC  205 

285  E\C  FILE  7 
END  FUE  7 

286  lv*l 
J*1 

ak«aki 

IF (KSXP) 536.531 .5^ 

531  fcRlTt  OUTPUT  Tap*  4.557 

557  FOWNAT  (1X34HT/*.  INPuLSF  RFSPCNSE) 


1 7 


290  HATIC(1>’»J)*T/4k 
G2*WAt!0<L^.J) *•? 

C 3*-G2/2. 

IF  (j5kP)  291 .297.201 

291  NP»Cn«l. 

If  IN)  ?96.29b,2')2 

292  00  29«.  L*  1 .IN 
NF»L*l 

NPRCC*NPRCC*NE.'.l 

2 95  CCMIM  E 

296  PRCDbM'ROU 

297  CFNaFRPD*PCP 
06*fc  XPF  ( Cj 3 ) 

CCEF2*G6/DEN 

G5SHAT!C<LW.J) 

IF (G5) 300. 520.120 
300  06*-G5 

Gc  TC  350 
320  G6*G5 
350  HSLF*  (1  ) *0. 

HSUFM2)  «!• 

IF  IN  ) 4 1 0 .6 10  ' 360 
360  CC  400  K*3.NP2 
AKV3*K-3 
KV1*<-1 
K^2* K-2 

hSU^  IK  ) «06*HSUN  ( K V l ) - AKN3»HSLM  (KM?) 

600  CONTINUE 
610  G7«HSUF'(NP2) 

G8«G7<M*2 

625  «rSP  ILF'»  J>  *COEF  2»G« 

If  I T m AX -T  ) 500.430.430 
630  JaJM 
T*T*P7 
GC  TC  290 
500  M«J-1 
N4*LW 

WRITE  OUTPUT  TAPE  4,55fl,((RATir(LW,I)  ,PESP(LF*»I)  »I*1.N3)  ,LF"N6,N6) 
558  FCPF‘AT(1X.F10»5,5X.F10#5) 

535  WRlTt, OUTPUT  Tapf  A , 5 16  . ( < PE  5P  I L"  . !)  . I « 1 ,N3 ) «LMbN<,  .N6 ) 

5U  FCRVAT  (F1C.5) 

IF (AK-AKNAX) 496.2 37, 237 
695  Aic«A<*F12 
Lv«LM* 1 

Ja  1 
TsTi 

GC  TC  290 
237  LVil 

EnO  FILE  8 
ENC  FILE  8 
N 5»0 
An»AK  1 
1 a 1 

502  Cll*n7/ A*. 

T.M 

J*  1 

ri2*Ar>«iF  IU11  > 

503  C5UF‘  (1  N.J1  ao. 


6C6  FCA  t'TsTAt  / r T 

F i *A“t<iA*UFLT  A*FRACT  / 2*0 
F?:f 1 »T Al 

F ■»  = !r  ? -F'F*  I 

F^scr^F  ( fe 3 / D F C-) 
fail  a'P(Tau*R) 

FC'-  = R ( I ) *L4 
TcTF* (T-TAU) /A* 

C«udEbAL'SF  ( TCI)  ) 

PlL^C-1-  APE-1 .0 
IF  (Plt'N)  *50*«5c  *2 79 
279  F9=C.n 

c-c  TC  AO  I 
«5C  E5«GBAPt/Dl2 
NF5aF6 
I A*NF.6*1 
IP*NE5*2 

E6«GPAPE-RATIC(LM»IA) 

F7«F6/Pll 

E8«1.0-E7 

R5PN»F7*Rt5P(L',«  IH)  ♦F0*RESP  (LM.  I A) 

Eq*FC^*RSPN*6lO 
801  GSLM  <(_►•«  J>  =C,Sl,*M  I N,J)*F9 

tau*tau*pio 

1*1*1 

E 10*Ar5F (TAU) 

IF  (TT«F10)520*505,505 
520  GSI.N  < I P * J ) *G5L‘*  ( LP  • J ) / AK 
IF  (TF'AX-T)  540,540,825 
525  T*T*fl7 
J*J*1 
TAt.aTAUl 

I F ( SENSE  S*ItCm  5)9997,9999 

9997  Cc  9998  L**1»15 
ICLMP(LK) **0 

9998  CONTINUE 
ICUNT ( 16) =-6 
ICCNP ( IT) **0 
IPU*P C 18) *N5 
CALL  OlJMP(  IDt'PP) 

9999  GC  TC  503 

540  If (AK-AKNAX) 543,843*860 
543  A«  * AK  *P  1 2 
N3«J-1 
N4*LW 

URlTk  OUTPUT  TAPE  6,559,  I (GSUM(Lm,I)  ,I«1«N3)  *L^»N4,N4) 
559  FCRV AT ( F 10.51 
N5»N5* J 

IF (AK-AKkAX) 278,278,860 
278  L V b L y » l 
TxTt 

TALslAl  1 
1*1 

f-C  TC  *02 
86c  E\ C * ! I E 6 
FN^  FT!  E 6 
536  ^t*M*l 
n?*M»? 


H 


I f < ISkP) 650.570,650 
570  CALL  Pt  0TS(RLFFER<102A)  .102A.5) 

CC  560  Jal.Nl 
XAXIS»J)S2<J> 

YAXjijj)  atfERAD(J) 

560  CONTINUE 

CALL  PLOT  (0.0.S.0.-5) 

call  scale  <yaxis.5.o.ni.i»io.O) 

CALL  SCALE  (XAxlS. 10. 0.N1. 1,10.0) 

CALL  LINE  UAXIS.YAXI5. Ml. 1,1*11) 

CALL  AX  1 5 (0.0.0.0.20H9EPACIATICN  FUNC TICN, 20*5.0. 90. 0»YAXIS(*1).Y 
1 AX  I S (Y2 ) .10.0) 

CALL  AXIS  (0.0,0.0.1H2.-1.10.0.0.0.XAXIS(y1) «XAXIS(M2) .10.0) 
call  plot  <0.0.0.0.9991 

650  READ  INPUT  T APf-  3.9502.E0 
IF  <ED-TEOC) 950  3.9509. 9 5C3 

9509  WRITt  OUTPUT  TAPE  A .95 1 1 
95ii  format (ahoenD) 

END  FILE  A 

9510  STOP  5 

END (1*1* 0.1*1) 


